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We modify the treatment of Multiple Parton Interactions (MPI) in Pythia by includ¬ 
ing the 1 0 2 mechanism and treating the 2 0 2 mechanism in a model-independent way. 
The 2 0 2 mechanism is calculated within the mean field approximation, and its parameters 
are expressed through Generalized Parton Distributions extracted from HERA data. The 
parameters related to the transverse parton distribution inside the proton are thus indepen¬ 
dent of the performed fit. The 1 0 2 mechanism is included along the lines of the recently 
developed formalism in perturbative QCD. A unified description of MPI at moderate and 
hard transverse momenta is obtained within a consistent framework, in good agreement 
with experimental data measured at 7 TeV. Predictions are also shown for the considered 
observables at 14 TeV. The corresponding code implementing the new MPI approach is made 
available. 
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I. INTRODUCTION 

It is widely realized now that hard Multiple Parton Interactions (MPI) play an important 
role in the description of inelastic proton-proton (pp) collisions at high center-of-mass energies. 
Starting from the eighties m until the last decade [5H26], extensive theoretical studies have been 
carried out. Attempts have been made to incorporate multi-parton collisions in Monte Carlo (MC) 
event generators [2MI]. Multiple parton interactions can serve as a probe for non-perturbative 
correlations between partons in the nucleon wave function and are crucial for determining the 
structure of the Underlying Event (UE) at Large Hadron Collider (LHC) energies. Moreover, 
they constitute an important background for new physics searches at the LHC. A large number of 
experimental measurements have been performed at the Tevatron |32H34j and at the LHC |35H38] . 
showing evidence for MPI at both soft and hard scales. This latter case is usually referred to as 
“Double Parton Scattering” (DPS), which involves two hard scatterings within the same hadronic 
collision. The cross section of such an event is generally expressed in terms of the (Teg. In the mean 
field approximation cje// [TI 1221 12^ , is the effective area which measures the transverse distribution 
of partons inside the colliding hadrons and their overlap in a collision. 

Recently, a new approach based in perturbative Quantum Chromodynamics (pQCD) has been 
developed HTHSi! for describing the MPI and its main ingredients are: 

• the MPI cross sections are expressed through new objects, namely double Generalized Parton 
distributions (GPD 2 ); 

• besides the conventional mean field parton model approach to MPI, represented by the so- 
called 2 (8) 2 mechanism (see Fig. 1 left), an additional 1 0 2 mechanism (Fig. 1 right) is 
included. In this mechanism, which can be described in pQCD, the parton from one of the 
nucleons splits at some hard scale and creates two hard partons that may participate in MPI. 
This mechanism leads to a significant transverse-scale dependence of MPI cross sections; 

• the contribution of the 202 mechanism to GPD 2 is calculated in a mean field approximation 
with model-independent parameters. 

The use of this new formalism at LHC experiments needs its implementation in MC event 
generators which has not been performed yet. The purpose of the present paper is to make a step 
ahead towards the implementation of this formalism into MC generators. We use the standard 
simulation of the MPI implemented in Pythia [5^, but with values of o'e// calculated by using 
the QCD-based approach of dM], i.e. including 10 2 processes. 
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FIG. 1: Sketch of the two considered MPI mechanisms: 2 0 2 (left) and 10 2 (right) mechanism. 


The current approach used for the description of the MPI in Pythia is based on |28l 129] . 
The Pythia code uses parton distribution functions, dependent on the impact parameter of the 
collision. From the theoretical point of view these are just one-particle Generalized Parton Distri¬ 
butions GPDi (see e.g. [401141] for a review). The parameters set in the Pythia simulation relative 
to the transverse parton density are extracted from fits to experimental data on UE, sensitive to 
the contribution of the MPI. This procedure is closely related to mean field based schemes, see 
e.g. [IT]. 

Such an approach has, however, a number of difficulties, both conceptual and practical. First of 
all, a problem arises at the level of mean field approximation. The transverse parton distributions 
have been extracted from J/T photoproduction measurements at the HERA collider, using QCD 
factorisation theorems da ng ioi E]. Hence they can not be treated as free parameters of the 
model. Secondly, it has been observed that different Pythia parameters are obtained when data 
sensitive to a different region of the MPI spectrum are used for the fits. For example, it has been 
shown [42] that different parameters result for fits to UE or hard MPI data. This might be an 
indication that an additional transverse scale dependence, which is not present in the mean field 
approach, is needed to describe experimental data on UE and hard MPI simultaneously. Recent 
improvements in the Pythia MPI model include a dependence of the parton transverse density 
on the longitudinal momentum fraction (x) |28|, but this only accounts for the x values of the 
hardest dijet. A complete x dependence which considers soft and hard partons may be irrelevant 
for the UE description where the transverse scales are rather close, but may become important for 
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measurements sensitive to hard MPI. 

The approach used in this paper combines the standard Pythia MPI model with the one of 
[TTHSni. We use a single gaussian to model the matter distribution function of the protons in 
Pythia. With these settings, the value of Cg// would be constant and independent on the scale. 
In order to implement the x and the scale dependence of Ue// in collisions where a hard MPI occur, 
these events are rescaled according to: 


^e// 


a 


(0) 

eff 


l + R’ 


( 1 ) 


where is the effective cross section in the mean field approach calculated in a model 

independent way from GPDi, parameterized from HERA data dMH, and R corresponds to 
the correction due to 1 (8> 2 mechanism [m i2oj. Such an approach is equivalent to using the 
GPDi-based transverse parton densities for parton transverse distributions. 


The main result of this paper is that the approach discussed above gives a unified description 
of both hard MPI and UE experimental data, with good accuracy and few fit parameters. The fit 
parameters are related to the amount of simulated MPI and of the color string reconnection, and 
to the separation scale between soft- and hard-scale processes, Qq whose value is expected to lie 
in the range 0.5-2 GeV^. The transverse scale dependent function R is calculated numerically by 
solving the nonlinear evolution equation |191 [20] . Predictions using this approach are shown later 
in the paper, and are labeled as “UE Tune Dynamic cr^ff”- Our analysis shows that the values 
of observable for UE are quite close to the results obtained in a free parton model (mean field 
approximation), while the inclusion of transverse scale dependent rescaling calculated in pQGD 
|10| improves the description of hard MPI. 

The paper is organized as follows. In section 2, basic theoretical ideas of the used approach are 
presented, while in section 3 their MG implementation is discussed. In section 4, comparisons for 
various predictions to observables measured at 7 TeV are shown. In section 5, predictions for these 
observables are presented for pp collisions at 14 TeV. In section 6 we compare our approach with 
the recently developed HERWIG-EE-5 approach, before drawing the conclusions in Section 7. 
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II. A SUMMARY OF THE THEORETICAL BACKGROUND 


The MPI four-jet cross section is characterized by the cross section cr^ff, which corresponds to 
an effective interaction area m, and can be written as: 


^^four-jet ^^two-jet ^^two-jet 


1 


( 2 ) 


dtl2dt3i dtl2 dt34 aeff 
where partons 1 and 2 create the first (12), and partons 3 and 4 the second (34) dijet. The pQCD 
calculation leads to the following expression for cJe// in terms of two-particle GPD: 


1 


^eff 


d?K 


[2]G2{xi,X3, Ql, Ql; A)[ 2 ]G 2 (x 2 , X4, Ql, Q\\ -A) 


+ [1]G2(xi, X3, Q\, QI; A)[2]G(x2, X4, Qi, Ql] -A) 

+ [i]G2{x2,X4,Ql,Ql]A)[2]G2ixi,X3,Qi,Ql]-A)]. (3) 


The second and third terms in Eq. correspond to the 1 (8> 2 mechanism, when two partons 
are generated from the splitting of a parton at a hard scale after evolution, while the first term 
corresponds to the conventional case of two partons evolving from a low scale, namely the 2 (8) 2 
mechanism. This first term can be calculated in the mean field approximation dSHIZ]. The 
momentum A is conjugated to the relative distance between the two participating partons. The 
fnll double GPD is a sum of two terms: 


G2{xi,X3,QI,QI,A) =[i] G2ixi,X3,Qj,Ql,A) -b[2] G2{xi,X3,QI,QI,A). (4) 

Here [ 2)^2 corresponds to the part of double GPD 2 , when both partons are evolved from the 
initial nonperturbative scale, while [i]G 2 corresponds to the case when one parton evolves up to 
some hard scale, where it splits to two successive hard partons, each of them in turn participating 
to the hard dijet event. We refer the reader to dZlIIE] for the detailed definitions of [ 1)^2 and 
[ 2 ]G 2 and their connection to light cone wave functions of the nucleon. 


For the two-parton GPD 2 we have: 

[2]G2{xi,X3, QI, QI, A) = D{xi, Qi)D{x3, Q 2 )F 2 g{A, xi)F 2 g{A, X3), 


( 5 ) 


where D{x, Q^) is a conventional parton distribution function (PDF). The use of the mean field 
approximation results in: 


l2]G2ixi,X3,Q‘l,Ql,A) = Gi{xi,Ql, A)Gi{xi, Ql, A) 


(6) 
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and 


Gi{xi,Ql,A) = D{xi,Qi)F2g{A,xi). 


(7) 


For the two gluon form factor F 2 g, we nse the exponential parametrization [16]. In fact, it leads 
to the same numerical results as the dipole form m, but it is more convenient for calculations. 
This parametrization is unambiguously fixed by J/T diffractive charmonium photo/electro pro¬ 
duction at HERA. The functions D are the conventional nucleon structnre functions and F 2 g can 
be parameterized as: 


F 2 g{A,x) = exp(-Rg(x)A^/2), 


( 8 ) 


where Bg{x)= Bq + 2^”^ dog(xo/x), with xq ~ 0.0012, Bq = 4.1 GeV“^ and Kq = 0.14 GeV“^. In 
onr implementation the central values of the parameters Bq and Kq [16] have been used, which are 
known with an accuracy of ~ 8 %. Integrating over A^, we obtain for the part of a^ff corresponding 
to the first term in Eq. 


^(0) 27r Bg{xi) + Bg{x2) + Bg{xQ) + Bg{xi)' 

J 

where X 1..4 are the longitudinal momentum fractions of the partons participating in the 2 0 2 
mechanism. This cross section corresponds to the free parton model and is model independent 
in the sense that its parameters are determined not from the fit of experimental LHG data, but 
from the fit of single parton GPDi. The maximum transversality kinematics i.e. 4Q^ = X 1 X 2 S for 
each dijet, have been considered in our approach, being Q the dijet transverse scale, and xi, X 2 the 
Bjorken fractions of the jets. 

The second and third terms in Eq. [^are parameterized as: 


^e// 


a 


(0) 

eff 


1 + R’ 


( 10 ) 


where R{Qi,Q 2 ,Qq) is calculated by solving iteratively the nonlinear evolntion eqnation, as ex¬ 
plained in detail in [I9l[20]. According to the results of m, the dependence of R on Xj in the 
maximum transversality regime is very weak and can be neglected with high accuracy. The function 
R also depends on the physical parameter Qq which corresponds to the separation scale between 
soft and hard dynamics where the GPD 2 is assumed to factorize. 





III. MONTE CARLO IMPLEMENTATION AND DEFINITION OF EXPERIMENTAL 

OBSERVABLES 
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In this paper we carry out two types of simulations: one based on the new approach defined in 
the sections 1 and 2 and one which follows the standard Pythia approach, used for comparison. 

Let us recall the standard Pythia approach which is referred as to ”UE tune” hereafter. In 
this study we use the Pythia 8.185 Monte Carlo event generator m It simulates a 2—7-2 matrix 
element interfaced to parton shower and Underlying Event (UE). The Pythia 8 event generator 
uses a simulation of the parton shower ordered in transverse momentum and the Lund string model 
|45j to implement the hadronization process. The performed study has considered as a starting 
point the UE simulation implemented in the Pythia 8 tune 4C |29j . This simulation makes use 
of the CTEQ6L1 [l8] parton distribution function and of a simple gaussian as a transverse matter 
distribution function. A fit to experimental data sensitive to the UE is performed in order to 
optimize the parameters related to the amount of MPI and colour reconnection in the simulation. 
The fit operation has been carried out by using the RIVET [l3] software, combined with the 
PROFESSOR machinery |44j . For the tune, two different observables have been considered at a 
center-of-mass energy of 7 TeV measured by the ATLAS experiment |39j . They are related to the 
multiplicity, ^chg, and the sum of the transverse momentum, PrmT, of the charged particles in the 
region transverse to the direction of the leading charged particle in each event. The performed fit 
has used only the data points corresponding to transverse momenta of the leading charged particle 
between 2.0 and 15.0 GeV. The exclusion of the very low pT region (< 2 GeV) is motivated by 
the fact that processes at those scales are expected to be dominated by soft physics, including 
diffractive processes and soft nonperturbative correlations, i.e. along the lines of m- The upper 
cut off is arbitrary, since its variation starting from 5 GeV does not change the values of the 
observables. 

The result of the fit consists of a new set of UE parameters implemented in the “UE tune” 
hereafter. The values of the Pythia 8 parameters obtained for the “UE tune” after the fit are 
shown in Table HI 

The first parameter listed in the table refers to the value of transverse momentum, p^, defined 
at ^/s = 7 TeV, used for the regularization of the cross section in the infrared limit, according to 
the formula 1/p^ —l/(pj- The second parameter is the probability of color reconnection 

among parton strings. The value of a^ff is found to be around 29.7 mb at 7 TeV; this value is 
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Pythia 8 Parameter 

Value obtained for the UE tune 

Multiplelnteractionsip^Ref 

2.659 

BeamRemnants:reconnect Range 

3.540 

Reduced x? 

0.647 

Teff (7 TeV) (mb) 

29.719 

Teff (14 TeV) (mb) 

32.235 


TABLE I: Pythia 8 parameters obtained after the fit to the UE observables. The value of pTORef is given 
at a reference energy of 7 TeV. The values of the reduced and of aeff at 7 and 14 TeV are also shown in 
the table. 


significantly smaller than the one obtained by tuning the correlation observables of the four-jet 
scenario [38], which is around 19-21 mb. Note that the value of 29.7 mb is quite close to the one 
determined in mean field approach [I71120]. 

After fitting the UE observables for the “UE tune” determination, the considered predictions are 
also tested against measurements sensitive to the hard spectrum of the MPI. Measurements of such 
type have been conducted by studying correlations between outgoing objects in a proton-proton 
collision, for instance in four-jet final states measured at 7 TeV by CMS |38j . In this scenario, 
two dijets have been selected at different transverse momentum; two jets are required to have pT 
larger than 50 GeV and they are classified as “hard-jet pair”, while the so-called “soft-jet pair” is 
composed by the two other jets selected with px greater than 20 GeV. Two correlation observables, 
AS* and A^fj^px that are sensitive to DPS, have been considered. They are, respectively, the 
azimuthal angle between the two dijet planes and the pT balance between the soft jets and are 
defined as follows: 


AS = arccos 


PT (pairI ) • PT {pair 2 ) \ 

\pT{'pairi)\ X \pT{pair2)\) ’ 


( 11 ) 


> rel 


PT = 


\p. 


jeh 




' jet2 I 


( 12 ) 


Ip7 N + Ipt 

where pairi {pair 2 ) is the hard (soft) jet pair and jeti {jet 2 ) is the leading (subleading) soft jet 


Let us now move to the new approach, based on the dynamical pQCD-based formalism, de¬ 
scribed in sections 1 and 2. The x and scale dependence of CTe// has been implemented by reweight¬ 
ing on an event-by-event basis the Monte Carlo simulation in presence of a hard and moderate 
MPI. The X dependence is given by Equation]^ where x\^2 are taken as the longitudinal momentum 
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fractions of the partons participating in the hardest scattering, while X 3^4 refer to the longitudinal 
momentum fractions of the partons participating in the hardest MPI. The scale dependence is ex¬ 
pressed by Equation!^ where R takes for Qi and Q 2 the scales of, respectively, the hard scattering 
and of the hardest MPI. Different values of Qq have been considered in the range between 0.5 and 
2 GeV2. 

We considered both the case of moderate MPI (i.e. MPI at scales of several GeV), relevant for 
UE, and hard MPI. 

Eor UE we treat separately the events where there is only one hard scattering, which are not 
rescaled, and the events with additional hard MPI. Eor the latter events two approaches were 
checked. Eirst, we rescaled these events according to Equationsand IT, taking as Qi and Q 2 for 
the R function the scales of the two hardest scatterings. As shown in Section IV, the influence of 
this rescaling is very small (less than 5%), with respect to the standard Pythia ”UE tune”. This 
may be connected both with the small values of R obtained for UE, and with the fact that the 
ladder splitting is roughly taken into account for such scales by the large value of the parameter 
p?. ~ 2 GeV. 

Subsequently, the second approach tried was to rescale only MPI events starting from the scale of 
order 4-5 GeV. When we rescale only the MPI starting from this (or a higher) scale, UE observables 
are not affected at all. At the same time, with this approach we avoid possible double counting 
effects, since at these scales the regularization formula in Pythia represents an ansatz for higher 
twist effects, including MPI. Thus, while using Pythia, we can neglect rescaling of MPI in UE, 
fitting instead. With the current accuracy, any of these two approaches can be used, leading to 
identical numerical results. This is in agreement with the approach documented in m. where it 
was argued that at scales relative to UE the values of (Tg// are close to the ones calculated in the 
mean field approximation. 

We consider now the case of hard MPI, specifically DPS. Two different processes may produce 
four jets in the final states. The first one is the so-called Single Parton Scattering (SPS) where 
the four jets are emitted through the same chain while the second one is DPS where the two hard 
interactions produce one dijet each. A different event topology is expected from these processes: 
if the four jets are produced through SPS, a high correlation between the objects of the final state 
is present and this is reflected in their relative configuration in the transverse plane. The direction 
of the hard jets, for example, is randomized by the emission of the additional two jets within the 
same chain and their initial pT balance is ruined. Instead, jet pairs coming from DPS events, 
namely from two independent scatterings, tend to be uncorrelated and their initial back-to-back 
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configuration is less subject to smearing effects coming from additional hard radiation; the jet pairs 
are expected to exhibit a more balanced configuration in pT and azimuthal angle. In particular, 
as shown in [SH], DPS events add a relevant contribution at low values of AS and Here 

we consider the experimentally relevant example, when the two dijet scales are 50 and 20 GeV. 

Similarly as before, the x and scale dependence of erg// have been implemented by reweighting 
on an event-by-event basis the Monte Carlo simulation, as explained above. 

In case that only MPI with px scales smaller than 15 GeV are present in the collision, no x and 
scale dependence is applied to the cJe// value of the corresponding event. The choice of 15 GeV is 
motivated by the fact that we need to treat differently the two contributing processes, SPS and 
DPS. Events where the two dijets are produced through SPS accompanied by moderate MPI, 
should not be reweighted piniiiMH]. In case a hard MPI occurs in the collision, dynamical agf / 
values are used. In this way, we assume that all collisions with a MPI scale greater than 15 GeV 
produce the second hard dijet px > 20 GeV pair selected in the considered four-jet scenario, while 
MPI at lower scales are below threshold for producing jets with px > 20 GeV. This approach is 
generally followed by standard experimental measurements for (Jeff determination, as the ones 
documented in [35l|36]. For our studies, lowering the 15 GeV cut off by 5-10 GeV shows variations 
of the predictions of DPS-sensitive observables of less than 2%. This is a clear indication of the 
consistency of our approach. 

Various simulation settings have been considered for comparison: 

• “UE tune”: predictions obtained with the parameters listed in Table ||] and without applying 
any reweighting of the simulation; this tune uses a constant value of cr^ff, following the 
standard Pythia approach; 

• “UE tune Q^-dep”: predictions obtained with the UE parameters listed in Table and by 
applying the scale dependence of a^ff with Qq = 1 GeV^; 

• “UE tune x-dep”; predictions obtained with the UE parameters listed in Table ||] and by 
applying the x dependence of cr^ff; 

• “UE tune Dynamic CTef/'' predictions obtained with the UE parameters listed in Table |I] 
and by applying both the x and the scale dependence with Qq = 1 GeV^. 

For the considered “UE tune Dynamic predictions using Qq values equal to 0.5, 1 and 2 

GeV^ have been also tested and compared. 
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A full MC implementation of the presented approach may be different from the one used in this 
paper, which relies on reweighted events simulated by Pythia. There are at least three reasons 
for it: 


• by using the Pythia event generator, all ladders are assumed to evolve independently from 
the low transverse scale, the initial-state radiation (ISR) being regularized by primordial 
gluon distribution with a transverse scale equal to No parton ladder splittings are 
included in this approach; 

• in Pythia, the geometric picture of the collisions in the impact parameter space corresponds 
to the 2(8)2 mechanism, while for 1(8)2 mechanism the geometrical picture would be different; 

• for multi MPI events, namely for events with several MPI within the same collision, we 
neglect the change of relative weight of 1 (8) 2 and 2 (8) 2 mechanisms. 

In this paper, these effects are neglected. First, the good agreement with experimental data 
shows that the high regularization scale may be a good alternative parametrization of the ladder 
splittings and of the corresponding changes in ISR at the UE transverse scales. In other words, 
for UE the high p^, which regularizes the charged particle multiplicity, also approximately fits the 
change of multiplicity due to ladder splitting. The Pythia regularization formula in this case can 
be viewed as an ansatz for twist expansion, that may include part of the MPI. Note that the ladder 
splitting scale is much smaller than the scales of hard dijets created by partons that evolve after 
the splitting [201 El] • So the effective ladder splitting is partly taken into account for UE by a high 
p^ value. This is the reason why the UE observables change only slightly in the new approach. On 
the other hand for hard MPI, when the hard splitting scale is much larger than p^, the inaccuracy 
in accounting for ISR at small pT can be safely neglected. 

Second, the direct calculation along the lines of [7| shows that neglecting the change of geomet¬ 
rical picture and of the relative weight between mean field and 1 (8) 2 mechanisms, when more than 
two separate dijet events are present, does not lead to numerical changes. 

We conclude that using events simulated with Pythia and reweighted with x- and scale- 
dependent values of (T^ff is a good approximation. In this way, we investigate the influence of 
changes of a^ff on MC observables sensitive to UE and DPS. 
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IV. RESULTS FOR 7 TEV 

In this Section, comparisons between UE- and DPS-sensitive measurements at 7 TeV and var¬ 
ious predictions are shown. Figure shows comparisons to ATLAS data [M| on charged particle 
multiplicity and px sum in the transverse region as a function of the leading charged particle px- 
Note that these are the observables which have been used in the fitting procedure for the determi¬ 
nation of the “UE tune”. The measurement is well reproduced by all considered predictions with 
discrepancies of only up to 10% in the high-p^ region (px > 10 GeV). The intermediate px region 
{2 < px < 10 GeV) is very well reproduced, while all predictions underestimate the first bins at 
Px > 2, GeV. This effect might be due to a not optimal simulation of diffraction in Pythia 8. 
However, no relevant differences are observed for the different a^ff models. 




——'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—' u.u '—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—'—' 

2 4 6 8 10 12 14 i6 i8 20 2 4 6 8 10 12 14 16 18 20 

p± (leading track) [GeV] px (leading track) [GeV] 


FIG. 2: Charged particle density (left) and prsum density (right) as a function of the leading charged particle 
in the transverse regions, measured by the ATLAS experiment at 7 TeV [33] ■ The data are compared to 
various predictions: the UE tune with constant CTe// value (red curve), the UE tune with Ue// x dependence 
applied (blue curve), the UE tune with a^ff scale dependence with Qq=1.0 GeV^ applied (black curve) and 
the UE tune with both cTeff x and scale dependence with Qq=1.0 GeV^ applied (pink curve). The lower 
panel shows the ratio between the various prediction and the experimental points. 


In Figure predictions obtained with different values of the scale Qq values are shown. All 
predictions are able to reproduce the measurement at the same good level. From this study, one 
may conclude that the UE data are not sensitive to the different settings of dynamical dependence 
applied to Ue//. 

Figure shows predictions with the varions Ue// settings considered previously, compared to 
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FIG. 3: Charged particle density (left) and prsum density (right) as a function of the leading charged particle 
in the transverse regions, measured by the ATLAS experiment at 7 TeV [33] ■ The data are compared to 
various predictions obtained with the “UE tune” where both the x and scale dependence have been applied 
for a^ff with Qo equal to 1.0 (red curve), 0.5 (blue curve) and 2.0 (black curve) GeV^. The lower panel 
shows the ratio between the various prediction and the experimental points. 


the normalized cross section distributions as a function of the correlation observables, AS and 
^ToftPT: measured in four-jet scenarios |38|. For these variables, the considered models show 
relevant differences. The static a^ff dependence (“UE tune”) is not able to properly describe the 
distribution as a function of AS; in particular, the region at low values (AS < 2.5), where a DPS 
contribution is expected, is underestimated by about 10-18%. By introducing the x dependence 
for cje// (“UE tune x-dep”), the agreement at low values of AS does not significantly improve. 
When the scale dependence of fie// is introduced (“UE tune Q^-dep”), the description of the 
normalized cross section as a function of AS gets better with differences not larger than 10%. The 
best agreement with the measurement is obtained for predictions where both the x and the 
dependence (“UE tune Dynamic aef/’) is included. The normalized cross section as a function of 
reproduced by all considered predictions. However, it has been already observed 
in jSS] that is less sensitive to a DPS contribution than AS, which uses information from 

both jet pairs. 

In Eigure predictions obtained with three different values of Qq (0.5, 1.0 and 2.0 GeV^) 
are compared to the normalized cross section distributions as a function of AS and A 

considerable level of agreement for the different settings is obtained. Predictions obtained with 
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FIG. 4: Normalized cross section distributions as a function of the correlation observables AS (left) and 
(right) measured in a four-jet scenario by the CMS experiment at 7 TeV |3S]. The data are compared 
to various predictions: the new UE tune (red curve), the new UE tune with the x dependence applied (blue 
curve), the new UE tune with only the scale dependence with Qq=1.0 GeV^ applied (black curve) and the 
new UE tune with both x and scale dependence with Qq=1.0 GeV^ applied (pink curve). The lower panel 
shows the ratio between the various prediction and the experimental points. 


Qq = 0.5 GeV^ are in good agreement with the AS measurement but overestimate the first bin of 
^^softPT- For Qq = 1 GeV^ and 2 GeV^ the agreement tends to improve for but is worse for 

AS. However, the measurement of the four-jet correlation observables is not able to discriminate 
the best choice for the value of Qq. 

In order to isolate the DPS contribution from the background produced by 2—)-4 processes, a 
dedicated event simulation has been performed with Pythia 8. Events with two hard scatterings 
within the same pp collision are simulated: the first hard scattering is generated with an exchanged 
transverse momentum between the outgoing partons, pT, larger than 45 GeV while for the second 
one, pt is required to be greater than 15 GeV. Figure shows the absolute cross sections predicted 
by the different settings implemented in the Pythia 8 simulation. 

The red curve shows the predictions for the UE tune with a static value of a^ff while the blue, 
black and pink lines represent the predictions obtained when implementing the dynamical x and 
dependence with y equal to 0.5, 1.0 and 2.0 GeV^, respectively. The highest DPS contribution 
is observed for Qq = 0.5 GeV^ and it decreases for increasing Qq values. The lowest contribution is 
observed for the static UE tune when no x and dependence is applied. The difference between 
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FIG. 5: Normalized cross section distributions as a function of the correlation observables AS (left) and 
(right) measured in a four-jet scenario by the CMS experiment at 7 TeV |3S]. The data are compared 
to various predictions obtained with the new UE tune where both x and scale dependence have been applied 
with Qg equal to 1.0 (red curve), 0.5 (blue curve) and 2.0 (black curve) GeV^. The lower panel shows the 
ratio between the various prediction and the experimental points. 
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FIG. 6: Absolute cross section distributions as a function of the correlation observables AS (left) and 
Ag^yjPT (right), produced via double parton scattering in a four-jet scenario at 7 TeV. Various predictions 
are shown in the figures: the new UE tune (red curve), the new UE tune with the x dependence applied 
(blue curve), the new UE tune with only the scale dependence with (5p=1.0 GeV^ applied (black curve) and 
the new UE tune with both x and scale dependence with Qg=1.0 GeV^ applied (pink curve). The lower 
panel shows the ratio between the various predictions and the predictions obtained with the new UE tune. 
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the static and the dynamical cJe// tune with Qq = 0.5 GeV^ is around 80%. The different DPS 
contributions observed among the considered predictions reflect the decreasing cTeff values for 
decreasing Qq as a function of the scale of the two scatterings (see Appendix of this paper). No 
significant differences in the shape of these distributions as a function of Qq are obtained. 

We observed that predictions of a dynamical cie// tune including a x and scale dependence of the 
transverse parton distribution are fully consistent with experimental data sensitive to moderate 
and hard MPI. The good agreement obtained for hard MPI is achieved due to contribution of 1 (8> 2 
mechanism. The contribution from this mechanism is essentially model independent, except for 
Qq [ 20 ] 5 which is the only new fit parameter, which is expected to lie in the 0.5 — 2 GeV^ range. 


V. PREDICTIONS FOR 14 TEV 


The dynamical a^ff dependence has been tested for predictions of UE and DPS observables 
at a center-of-mass energy of 14 TeV. The x and scale dependence of Ue// follows respectively 
Equations and 10 similarly as for 7 TeV. Note that the function R in Equation |10| also depends 
on the center-of-mass energy y/s |TH|. Figure shows predictions of charged particle density and 
the pt sum as a function of the leading charged particle pT, while in Figure the normalized cross 
sections as a function of the four-jet correlation observables, AS and A^fj^^pT, are presented. The 
predictions have been obtained by using the UE tune with a static cJe// value and with a dynamical 
X- and Q^-dependent a^ff value, with various values for Qq: 0.5, 1.0 and 2.0 GeV^. 

For each plot, the ratio to the predictions obtained with the UE tune with a constant cr^ff 
value is shown. While for the considered UE observables, a very small change is observed for 
the various predictions, larger differences are observed when the four-jet correlation observables 
are investigated. In particular, tunes with a dynamical Ue// dependence tend to predict a higher 
contribution at low AS and A^^^^ppT values. These are the regions where a contribution from DPS 
is expected. The difference between static and dynamical dependence is of up to 15% for 
AS < 2.0. Predictions with Qq equal to 1.0 and 2.0 GeV^ are very similar to each other, while 
results obtained with Qg = 0.5 GeV^ show a higher contribution at low values of AS and 
where the contribution of hard MPI is expected to be relevant. 
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FIG. 7: Charged particle density (left) and prsum density (right) as a function of the leading charged 
particle in the transverse regions at 14 TeV. Various predictions are shown in the figures: the new UE tune 
(red curve), the new UE tune with both x- and scale-dependence with Qq=1.0 GeV^ (blue curve), Qq=0.5 
GeV^ (black curve) and Qq=2.0 GeV^ applied (pink curve). The lower panel shows the ratio between the 
various prediction and the experimental points. 
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FIG. 8: Normalized cross section distributions as a function of the correlation observables AS (left) and 
(right) in a four-jet scenario at 14 TeV. Various predictions are shown in the figures: the new 
UE tune (red curve), the new UE tune with both x and scale dependence with Qq=1.0 GeV^ (blue curve), 
Qq= 0.5 GeV^ (black curve) and Ql=2.0 GeV^ applied (pink curve). The lower panel shows the ratio between 
the various predictions and the predictions obtained with the new UE tune. 
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VI. COMPARISON WITH RECENT HERWIG TUNES 

The calculations described so far in this paper are based on the MPI approach implemented in 
Pythia. a different approach for the description of MPI is implemented in the Herwig++ event 
generator |23H25j . Recently, a new tune has been released for the simnlation of the UE, labelled 
as UE-EE-5-CTEQ6L1 |24j . This tune is very interesting for the purpose of this paper because it 
is able to simultaneously describe data sensitive to soft MPI and predict a value of Ue// of about 
15 mb, which is much lower than the one in ’’Pythia UE tune”. However, the approach of the 
UE-EE-5-CTEQ6L1 tune is based on a very different picture of both UE and hard MPI than the 
one discussed in our paper: 

• in [231I25| . the mean field approximation is used to describe hard MPI, with parameters 
related to the transverse parton density distribution obtained through a fit to hard MPI 
data. The parametrization of the transverse parton distribntion corresponds to a dipole 
form of the two gluon form factor (Eq. equal to: 

^ (l + AVm^) • 

The parameter |231l25j has the same physical interpretation as the parameter in¬ 
troduced in [niiiT!, measuring the gluonic radius of the proton. In “UE Tune Dynamic 
(Te//” developed in this paper, the transverse parton distributions have been determined 
from HERA data [T5HT7] . having thus the parameter as a model-independent input. 
Comparing and m^, i.e. comparing the values of the gluonic radii used by tunes UE-EE- 
5-CTEQ6L1 and “UE Tune Dynamic cTe//”, respectively, one gets ~ 2'm?g. This means 
that in the UE-EE-5-CTEQ6L1 tnne the gluonic radius of the proton in hard MPI is ^/2 
times smaller than the one observed in HERA. In our approach the gluonic radins of the 
proton is compatible with the one observed at HERA, but in addition to the mean field 
approximation, a 1 (8) 2 mechanism is included. The contribution of the 1 (8> 2 mechanism to 
hard MPI is of the same order as of the mean field approximation; 

• in order to describe UE data and to predict Ue// around 15 mb, the UE-EE-5-CTEQ6L1 
tune nses a color reconnecgtion model developed in [25] . In such approach one gets the value 
of ~3.9 GeV for the regularization threshold, of the partonic cross section. Eor “UE Tune 
Dynamic Ug//”, the description of UE data and the corresponding parameters are similar to 
“Pythia UE tune”. In particnlar, the value of implemented in “UE Tune Dynamic CTeff” 
is ~2.68 GeV (see Table |^; 
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• the MPI model implemented in Herwig++ does not lead to any transverse dependence for 
the value of cr^ff, which is taken as a constant as a function of the scale of the secondary 
hard scattering, in difference to the current approach. 


Predictions of the two described Herwig++ tunes have been compared to data sensitive to 
hard MPI. Figure]^ shows predictions of the old UE-EE-4-CTEQ6L1 [2l| and UE-EE-5-CTEQ6L1 
tunes, compared to the normalized distributions as a function of the correlation observables, AS 
and measured by CMS in four-jet final states at 7 TeV jSH]. Predictions from both tunes 

do not give a good description of the experimental data; UE-EE-5-CTEQ6L1 tune performs better 
than UE-EE-4-CTEQ6L1 but differences of around 20-30% with the data are observed for values 
of AS smaller than 2.5. 


Normalized AS in pp^ 4j in \r]\ < 4.7 at ^/s = 7 TeV 
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FIG. 9: Normalized cross section distributions as a function of the correlation observables AS (left) and 
(right) measured in a four-jet scenario by the CMS experiment at 7 TeV [35] • The data are compared 
to predictions obtained with Herwig tune UE-EE-4-CTEQ6L1 and tune UE-EE-5-CTEQ6L1. The 
lower panel shows the ratio between the various prediction and the experimental points. 


In conclusion, the approach used by the “UE tune Dynamic CTe//” developed in this paper and 
by the Herwig-|— |- UE-EE-5-CTEQ6L1 tune [21], are rather different and are based on a different 
picture of both UE and hard MPI. In “UE tune Dynamic cr^ff”, the emerging treatment of UE is 
quite close to mean field approach based on transverse parton densities determined from HERA, 
and ladder splittings (1 (8> 2 mechanisms) become important in the description of processes with 
hard MPI. In the approach of Herwig-|— |- UE-EE-5-CTEQ6L1 tune, soft and hard MPI are both 
described in mean field approach, but with a gluon radius of about 1.4 times smaller than the one 
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obtained from exclusive diffraction measurements at HERA, and a new color reconnection model. 
No 1 (8) 2 mechanism is included. We believe, that additional experimental data sensitive to soft 
and hard MPI will be able in the future to further constrain and eventually discriminate the two 
approaches. 


VII. CONCLUSIONS 

We have developed a new tune “UE tune Dynamic cTe//” [50]. The code modifies the treatment of 
hard Multiple Parton Interactions (MPI) in Pythia 8, leading to an improvement in the description 
of experimental data. We do not change the Monte Carlo code of Pythia, but we rather use the 
results of the MPI simulation on an event-to-event basis, so that 1 (8> 2 mechanisms are included. 
The tune uses a fit to Underlying Event (UE) data in order to extract the parameters relative to soft 
MPI and includes values of which contain the 1 (8> 2 mechanism. They are calculated directly 
in the mean field + pQCD approach, as discussed in [20]. The dynamical dependence of a^ff is 
not derived from a process-dependent fit of the experimental data, but is directly obtained from 
theoretical calculations dZHSO]. For the parameter Qq, that separates soft and hard scales, we have 
considered a range of values 0.5 < Qq < 2 GeV^. At present, the accuracy of the experimental data 
does not allow to carry a more precise determination, although the central values of the measured 
observables are better described by 0.5 < Qg < 1 GeV^. We observe that predictions from such 
tune are in good agreement with experimental measurements at 7 TeV, and for the first time give 
a consistent description of MPI at both moderate (UE) and hard scales. The results for UE are 
close to mean field approximation values, as anticipated in [19]. The additional transverse scale 
dependence of erg//, relative to mean held approach, due to 1 (8 2 mechanism, is essential for a 
unihed description of UE and hard MPI. 

Predictions, obtained with the new tune for proton-proton collisions at 14 TeV, which are 
expected to happen within the next LHC phase, are also presented. 
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Appendix A: aeff dependence at different energies for various scale and longitudinal 

momentum fraction choices 
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In this Section, a closer look at the a^ff dependence on scale, longitudinal momentum fraction 


and collision energy is provided. Figure 10 shows the values of a^ff as a function of the scale of 
the 2’^'^ interaction for a scale of the first interaction equal to 50 GeV and different choices of Qq 
( 0.5, 1.0 and 2.0 GeV^). In this study, the longitudinal momentum fractions of the first interaction 
system has been set to 0.014, corresponding to the maximal transversality regime. The x value 
relative to the second hard scattering has been also fixed to the maximal transverse momentum 
exchange. One can see that a^ff spans over a range of values between 16 and 30 mb, depending on 
the choice of Qq. The value of cr^ff decreases as a function of the scale of the 2”'^ hard interaction, 
Q 2 , showing a difference of about a factor of 1.1-1.2 between Q 2 = 15 GeV and Q 2 = 40 GeV. A 
significant dependence of aeff on the choice of Qq is also observed. The smallest agff values are 
obtained for Qq=0.5 GeV^, while they increase of roughly a factor of 1.25 and 1.44, for respectively 
Qq=1.0 and Qq=2.0 GeV^. 



FIG. 10: Values of Ue// as a function of the scale of the 2"^* interaction for different scales of the first 
interaction, Qi, and different choices of Qq. The values of the longitudinal momentum fractions correspond 
to the maximal transverse momentum exchange. 


In Figure 11, the aeff dependence is studied for various scales of the first interaction (50, 100 






22 



FIG. 11: Values of tTe// as a function of the scale of the 2"'^ interaction for different scales of the first 
interaction, Qi. The value of Qq has been kept fixed to 1.0 GeV^. 



FIG. 12: Values of a^ff as a function of the scale of the 2”^^ interaction at different collision energies at 
7 TeV and 14 TeV for first hard interactions occurring at a scale Qi = 50 GeV. The three values of Qq 
equal to 0.5, 1.0 and 2.0 GeV^ are considered and the longitudinal momentum fractions of the two dijets 
correspond to the maximal transverse momentum exchange for both y/s = 7 TeV and y/s = 14 TeV. 
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and 200 GeV) corresponding to choices of xi and X 2 in the maximum transversality regime, equal 
to respectively 0.014, 0.028 and 0.056. The values of X 3 and 2:4 related to the partons participating 
in the secondary hard scattering are also set to the maximal exchanged transverse momentum. In 
this study, only predictions obtained with Qq = 1.0 GeV^ are shown. It is observed that Ue// does 
not show a large dependence on the scale of the first interaction: in particular, (Te// decreases as 
a function of the scale of the first hard scattering. The three curves are very similar between each 
other as a function of the scale of the second hard interaction and the difference is less than 1 mb. 

Figure [T^ considers the aeff variation at different collision energies, 7 and 14 TeV, as a function 
of the scale of the second hard interaction. The three values of Qq equal to 0.5, 1.0 and 2.0 GeV^ 
are considered. Only scales of the first interaction equal to 50 GeV are examined. The value of 
aeff increases for increasing collision energies. For Qq=0.5 and 1.0 GeV^, aeff increases of about 
2-3 mb for any scale of the second hard scattering, while for Qq = 2.0 GeV^, the increase of aeff 
is larger and it reaches values of up to 4.5 mb at Q 2 = 15 GeV. 
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